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Abstract 

The Chebyshev spectral collocation method for the Euler gas-dynamic 
equations is described. It is used with shock fitting to compute several two- 
dimensional gas-dynamic flows. Examples include a shock/acoustic wave 
interaction, a shock/vortex interaction, and the classical blunt body 
problem. With shock fitting, the spectral method has a clear advantage over 
second order finite differences in that equivalent accuracy can be obtained 
with far fewer grid points. 
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Nomenclature 


a 

b 

(k^,ky) 

k 

P 

(r,0) 
s( t) 


u 

u 

Un 

(u,v) 

(x,y) 

(xQ,yo) 

y£ 

A 

B, C 
L 

M, N 

P 

Q 


^pq 

q( 1 , 0) 

^pjq 


sound speed 

vortex softening length 

constants in discrete Chebyshev transform 

wavevector 

wavevector magnitude 

pressure 

physical polar coordinates 
blunt body boundary 
shock front radius 
start-up function for linear waves 
physical time 

start-up time for linear waves 
solution to 1-D test problems 
interpolating polynomial 
solution at collocation points 
discrete Chebyshev coefficient 
physical velocities 
collocation points 
computational left boundary 
shock front location 
physical Cartesian coordinates 
center of downstream vortex 
periodicity length in y 
pressure wave amplitude 
coefficient matrices in Euler equations 
discrete spatial operator 
number of collocation points 
shock Mach number 
logarithm of pressure 
vector of dependent variables in the 
Euler equations 
spectral coefficients of Q 

spectral coefficients of 
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Q 


( 0 , 1 ) 

pq 


spectral coefficients of 


Q 


y 


R 

R'^, R" 

s 


T 

(U,V) 

(X,Y) 

3 

Y 

Ax 

At 

1C 

p 

a 

0 ^ 


0 


1 


max 


source term In the Euler equations 
linearized 1-D Riemann variables 
entropy (divided by specific heat at 
constant volume) 
computational time 
contravariant velocity components 
computational coordinates 
vertical stretching parameter 
ratio of specific heats 
mesh size 
time increment 
vortex circulation 
density 

smoothing function 

filter cut-off angle 

incident angle of pressure wave 

angle of computational outflow boundary 

Chebyshev polynomial of degree n 

s treamfunct ion 


iv 



I* Introduction 


No convincing case has yet emerged for the spectral shock-capturing 
technique for the Euler equations. Although solutions can be obtained by such 
a method, their accuracy tends to be quite low. The filtering procedures 
necessary to control the oscillations arising from the discontinuity in the 
solution at the shock have the side-effect of reducing the accuracy in the 
structured regions of the flow. In a previous paper^ we reported our 
experience with a spectral shock-capturing method on a periodic, 
one-dimensional, compressible flow problem. We found the method to be only 
first-order accurate and certainly no better than finite difference solutions. 

In the present paper we propose a straightforward cure for the 
oscillations that plague spectral shock-capturing methods: resort to spectral 

shock-fitting methods instead. 


II. Spectral Methods for Shock-Fitting 

Shock-fitting techniques have been a standard finite difference tool for 
some 15 years or so. They are suitable for problems in which the general 
features of the solution (but not the details) are predictable. This approach 
overcomes the difficulties of shock-induced oscillations for the simple reason 
that the shock front itself is a computational boundary. While it eliminates 
the need for computing derivatives across the shock (which is the source of 
the oscillations), it adds the complexity of requiring an algorithm to 
determine the shape and motion of the shock. Spectral methods for shock- 
fitting are a straightforward combination of both standard techniques. They 
do, however, require the use of Chebyshev polynomials rather than Fourier 
series in at least one coordinate direction. We begin this section by 
discussing the fundamentals of Chebyshev spectral methods. 
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Basic Chebyshev Spectral Concepts 


Consider the model problem 


Ut + 


= 0 


( 1 ) 


on -1 < X < 1 with initial condition 


u(x,0) = sin(2.5 ttx) 


( 2 ) 


and boundary condition 


u(-l,t) = sin[2.5 ‘ir(-l-t)). 


(3) 


The expansion functions are the Chebyshev polynomials 


T (x) = cos(n cos ^ x) 
n 


(4) 


and the collocation points are 


X. = cos —*7 
J N 




(5) 


Note that 


T 


TTjn 


„(*j) =■ c°s H 


( 6 ) 


The discrete Chebyshev coefficients are 


2 

Ncn 


N 

I 


j=0 
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(7) 
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where 


2 n = 0 or N 

1 1 < n < N-1 

Thus the interpolating function is 

N 

/s 

u(x) = } U T (x) • 

n n 

n=0 

The analytic derivative of this function is 



9u 

3x 


N 

I 


n=0 


u 


( 1 ) 

n 


Tn(x), 


( 8 ) 


(9) 


( 10 ) 


where 


-^(1) 

u =0 
N+1 


( 1 ) 


= 0 


( 11 ) 


.(1) -^(1) 


u_ 
n n 


Un+2 + 2(n+l)u^^.l 


n=N-l,N-2,«**,0 


The Chebyshev spectral derivatives at the collocation points are 


9u 


= I “n 

• ^ n 


( 1 ) 


n=0 


cos 


( 12 ) 


Special versions of the Fast Fourier Transform (FFT) may be used for 
evaluating the sums in Eqs. (7) and (12). The total cost for a Chebyshev 
spectral derivative is thus 0(N £n N) . 
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The time-stepping scheme for Eq. (1) must use the boundary conditions to 
update Ujj (at x = -1) and the approximate derivatives from Eq. (12) to 
update Uj for j=0 , 1 , • • • , N-1 . Note that no special formula is required for 
the derivative at j = 0 (or x = +1). Results at t = 1 for a Chebyshev 
spectral method, a Fourier spectral method and a second-order finite 
difference method are given in Table I. (The temporal discretization errors 
are negligible in all cases.) For this non-periodic problem Fourier spectral 
methods are quite inappropriate, but the Chebyshev spectral method is far 
superior to the finite difference method. 


Table I. Maximum Error for a 1-D Dirichlet Problem 


N 

Chebyshev 

Spectral 

Fourier 

Spectral 

Finite 

Difference 

4 

1.49 

(0) 

1.85 (0) 

1.64 

(0) 

8 

6.92 

(-1) 

1.92 (0) 

1 .73 

(0) 

16 

1.50 

(-4) 

2.27 (0) 

1.23 

(0) 

32 

3.45 

(-11) 

2.28 (0) 

3.34 

(-1) 

64 

9.55 

(-11) 

2.27 (0) 

8.44 

(-2) 


The Chebyshev collocation points are the extreme points of x^(x) . Note 
that they are not evenly distributed in x, but rather are clustered near the 
endpoints. The smallest mesh size scales as l/N^. While this distribution 
contributes to the quality of the Chebyshev approximation and permits the use 
of the FFT in evaluating the series, it also places a severe time-step 
limitation on explicit methods for evolution equations. 
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Filtering for Chebyshev Spectral Methods 

The same types of filtering operations that were discussed in Ref. 1 for 
Fourier spectral methods are applicable to Chebyshev spectral methods as well. 
In the latter case, however, there is as yet no theoretical support for the 
usefulness of pre-processing or derivative filtering on simple linear 
problems. 

A straightforward filtering procedure is to mimic Eq. (8) of Ref. 1 by 
setting 


n=0 ^ ^ 


(13) 


where u^ is given by Eq. (7) and a(0) is a standard smoothing function as 
described in Ref. 1. There are two problems with this approach: boundary 
conditions and conservation properties. Neither survives under this type of 
filtering. The lack of conservation in filtering does not appear to be 
crucial. After all, the shocks are not being captured, and as will be evident 
below, the computations use the Euler equations in nonconservation form. Any 
drift in the mean flow properties in the calculations has been minor. The 
boundary conditions are another matter. They are enforced after every 
application of filtering. 


A Spectral Shock-Fitted Method 

A schematic of the type of spectral shock-fitted calculations described 
below is illustrated in Fig. 1. At time t = 0 an infinite, normal shock 
at X = 0 separates a rapidly moving, uniform fluid on the left from the 
fluid on the right which is in a quiescent state except for some specified 
fluctuation. The initial conditions are chosen so that in the absence of any 
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fluctuation the shock moves uniformly in the positive x-direction with a Mach 
number (relative to the fluid on the right) denoted by Mg. In the presence 
of fluctuations the shock front will develop ripples. The shape of the shock 
is described by the function Xg(y,t). The numerical calculations are used to 
determine the state of the fluid in the region between the shock front and 
some suitable left boundary x^Ct) and also to determine the motion and shape 
of the shock front itself. 

Figure 1 is taken from a shock/turbulence calculation^ in which the 

downstream fluctuation is a plane vorticity wave that is periodic in y with 
period y Because of the initial value nature of the calculation, the fluid 
motion behind the shock is not periodic in x, as Fig. 1 makes abundantly 

clear. The interesting physical domain is given by 

x ( t) < x < X (y, t) (14a) 

Li S 

0 < y < (14b) 

t > 0. (14c) 


The change of variables 


X = 


X - Xj^( t) 

Xg(y,t) - x^(t) 


(15a) 


Y = yly^ 


(15b) 


T = t 


(15c) 
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produces the computational domain 

0 < X < 1 

0 < Y < 1 (16) 

T > 0. 

The fluid motion is modeled by the two-dimensional Euler equations. In 
terras of the computational coordinates these are 

Qt + B Qx + C Qy " 

/ 

where Q = (P,u,v,S)^, 

U yx^ yx 

g2 

— X U 0 

y X 

g2 

— X 0 U 

Y y 

0 0 0 





The contravariant velocity components are given by 
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U = + uX + vX 

t X y 

and (20) 

V = + uY + vY . 

t X y 

A subscript denotes partial differentiation with respect to the indicated 
variable* Reference conditions at downstream infinity are used to normalize 
p and S; u and v are velocity components in the x and y directions, 
both scaled by the characteristic velocity defined as the square root of the 
pressure-density ratio at downstream infinity* A value y = 1 *4 has been 
used* 

Let n denote the time level and At the time increment* The time 
discretization of Eq* (17) is 

Q = [l - AtL^jq’^ (21) 

j [q" + (1 - AtL)q], (22) 

where L denotes the spatial discretization of B 8 + C 8 * The solution 

X I 

Q has the Chebyshev - Fourier series expansion 

M N/2-1 „ 

Q(X,Y,T) ^ I I Q„„(T) T (23) 

p=0 q=-N/2 P‘1 P 

where ? = 2X-1. The derivatives and Qy are approximated by 


Qv = 2 


M 

I 

p=0 


N/2-1 

I 

q =-N/2 


q(1,0)(t)^ (^)e2iTiqY 

pq p 


(24) 
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Q 


Y 


M N/2-1 
= 21T I I 

p=0 q=-N/2 


q(0,1) 

pq 


(T)t (5)e 
p 


2TriqY 


(25) 


( 1 , 0 ) 

where Qp^ 


is computed from 


<pq 


in a manner analogous to Eq« (11) and 


q(0,1) 

pq 


i q Q 


pq* 


(26) 


The most critical part of the calculation is the treatment of the shock 
front* The shock-f itting approach used here is desirable because it avoids 
the severe post-shock oscillations that plague shock-capturing methods. The 
time derivative of the Rankine-Hugoniot relations provides an equation for the 
shock acceleration. This equation is integrated to update the shock position 
(see Ref. 3 for details). This method is a generalization of the finite 
difference method developed by Pao and Salas^ for their study of the 
shock/vortex interaction. 


Boundary Conditions 

The correct boundary conditions at the left boundary depends upon the 
relative Mach number. For uniform flow and y “ 1 the flow behind the 

shock is supersonic if > 2.08. In this case the boundary at is a 

supersonic inflow boundary and it is appropriate to specify all variables. If 
< 2.08 then the left boundary is a subsonic inflow boundary. The 
advisable procedure, then, is to base the numerical boundary conditions on the 
linearized characteristics of the Euler equations. At the left (subsonic) 
boundary the (linearized) characteristic variables corresponding to the 
outgoing characteristic direction are 
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R = P - 1 u. (27) 

a 

Similarly, 

= P + - u (28) 

a 

corresponds to the outgoing characteristic direction at the right (subsonic) 
boundary which is used by the shock fitting algorithm, 

A set of successful boundary conditions on the left is obtained by first 
calculating preliminary values of all quantities at the left boundary and then 
incorporating the given values of S, v, and R**" as 


S = S . 

given 


V = V , 

given 


P +Iu 

a 


+ 

R . 

given 


P -lu 


= p 


a prelim 


-- 1 - 

a prelim 


(29) 


Thus, the PDE is used to update the appropriate characteristic combination of 
variables at the boundary. The characteristic analysis is given in Ref. 5. 
The particular boundary condition was advocated in Ref. 6. For the right 
boundary a similar characteristic correction procedure can be incorporated 
into the evaluation of the shock velocity. 

The global nature of spectral methods makes them even more sensitive to 
the boundary conditions than finite difference methods. An illustration of 
just how unforgiving spectral methods can be is provided in Fig. 2. Shown 
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there are two spectral shock-fitted calculations of the interaction of a Mach 
1.3 shock with a Karman vortex street. (See Ref. 7 for more details about 
this problem.) The top row shows what happens when all flow variables are 
specified at the left, subsonic, inflow boundary, whereas the bottom row 
displays a calculation which is identical except for the use of Eq. (29) as 
the inflow boundary condition. The former calculation is clearly contaminated 
by oscillations emanating from the inflow boundary. The latter calculation 
makes clear that no physical signals have yet reached the inflow location even 
though in the spectral method numerical signals reach the inflow 
instantaneously. Finite difference calculations for this same problem were 
reported in Ref. 7. Despite the fact that an overspecified inflow boundary 
condition was used, no analogous problem arose because of the local nature of 
the discretization. 


III. Results for Chebyshev Spectral Shock-Fitting 

Shock/Turbulence Interaction 

The nonlinear interaction of plane waves with shocks was examined at 
length in Ref. 2. The numerical method used there was similar to the one 
described above but employed second-order finite differences in place of the 
present Chebyshev-Fourier spectral discretization. Detailed comparisons were 
made in Ref. 2 with the predictions of linear theory.^ The linear results 
turned out to be surprisingly robust, remaining valid at very low (but still 
supersonic) Mach numbers and at very high incident wave amplitudes. The only 
substantial disagreement occurred for incident waves whose wave fronts were 
nearly perpendicular to the shock front. This type of shock-turbulence 
interaction is a useful test of the spectral technique because the method can 
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be calibrated in the regions for which linear theory has been shown to be 
valid. 

The most reliable numerical results can be obtained for the acoustic 
responses to acoustic waves. Unlike the vorticity responses, these require no 
differentiation of the flow variables, thus eliminating one extra source of 
error. Moreover, the acoustic reponse stretches much further behind the shock 
than the vorticity response, thus providing greater statistical reliability. 
Vorticity response results are reported in Ref, 9, The incident pressure wave 
is taken to be 

i( k, - CO, t) 

Pj = AJ e (30) 

where kj^ = “l~ ^s®l ^1 ''^l A^ is the amplitude. In 

terms of the incidence angle 0^ , kj = (k^ cos 9^ ,k^ sinSj). The linearized 
transmitted acoustic wave can be expressed in the same manner with all 
subscripts changed from 1 to 2. The amplification coefficient for the 

transmitted acoustic wave is then the ratio A' /A' . Figure 3 indicates the 
transmission coefficient extracted from the computation. At each fixed value 
of X we perform a Fourier analysis in Y of the pressure. The Fourier 
coefficient for q = 1 provides the amplitude A'. In order to reduce the 

transients that would accompany an abrupt start of the calculation at full 
wave amplitude, an extra factor of s(t) is inserted into Eq. (30), where 

(3(t/t )^ - 2(t/t )^ 0 < t < t 

s(t) = { ® ® ® . (31) 

1 t > t 

V S 

The start-up time tg is some multiple (typically V 2 ) the time it takes 
the shock to encounter one full wavelength (in the x-d irect ion) of the 
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incident wave. The ratio ^£^^1 plotted in Fig. 3 as a function of the 

mean value of the physical coordinate x corresponding to X. The start-up 
time for this Mach 3 case is tg = 0.56. The average of the x-dependent 

responses between the start-up interval and the shock produces the computed 
transmission coefficient. The standard deviation of the individual responses 
serves as an error estimate. 

The dependence upon incidence angle of the acoustic transmission 
coefficient for A' = 0.001, and Mg = 3 waves is displayed in Fig. 4. As is 
discussed in Ref. 2, linear theory is quite reliable at angles below, say, 
45^. Figure 4 contains results from both spectral and finite difference 

calculations. The finite difference results were obtained with the same 
second-order MacCormack's method that was described in Ref. 2 except that 
periodic boundary conditions (rather than stretching) were employed in the y- 
direction. The finite difference grid was 64 x 16 and these calculations used 
a CFL number of 0.70. The spectral grid was 32 x 8, the CFL number was 

0.50. (No solution smoothing was applied.) Figure 4 shows that both methods 
produce the same results. A head-to-head comparison of both methods for the 
= 10^ case is provided in Table II. The ”exact" value is taken from 
linear theory.® Since the amplitude of the incident acoustic wave is so 
small, it should come as no surprise that four points in the y-direction 

suffice for the spectral calculation. Note that the standard deviations are 
substantially smaller for the spectral method. These results suggest that the 
spectral method requires only half as many grid points in each coordinate 


direction 


Table II. Grid Dependence of Acoustic Transmission Coefficient 


Grid 


Finite 

Difference 

Chebyshev- 
Fourier Spectral 

16 X 

4 

6.403 

± 

2.652 

7.257 

± 

0.587 

16 X 

8 

6.427 

± 

2.626 

7.257 


0.587 

32 X 

4 

7.105 

± 

0.453 

7.158 

± 

0.022 

32 X 

8 

7.134 

± 

0.471 

7.158 

± 

0.022 

32 X 

16 

7.139 

± 

0.497 

7.158 

± 

0.022 

64 X 

16 

7.163 

± 

0.078 

7.157 

± 

0.017 

X 

00 

CM 

16 

7.152 

± 

0.022 




**exact 

fl 

7.156 



7.156 




Shock /Vortex Interaction 

This problem is closely related to the previous one* The downstream 
field is not the linear plane pressure wave of Eq. (30) but an idealized 
vortex in which the density is constant, the velocities are derivable from the 
stream function 


^ log /b^ + (x - + (y - y^)^ , (32) 

the pressure from Bernoulli's equation, and the temperature from the equation 
of state* This model approaches an idealized incompressible point vortex at 
large distances from the vortex center at (xQ,yQ), but it is much smoother in 
the core* The specific example provided here has the circulation k = 0*40; 
the vortex softening scale b = 0*1 and the vortex is located at 
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Note that periodic boundary conditions in y are no longer appropriate* 
Accordingly, Eqs. (14b) and (15b) are replaced by 

—00 ^ y ^ 

and 

Y = tanh(3y) + 1 
2 

respectively, where the stretching parameter 3 is of order one* Moreover, 
the spectral method now uses Chebyshev series in Y as well as X* The 
analogs of Eqs* (14) - (26) are given in full in Ref* 7* (incidentally, it 
was this Chebyshev - Chebyshev algorithm which was used in the production of 
Fig* 2*) 

The computed results for the shock-vortex interaction at t = 0*35 are 
given in Fig* 5 for both finite difference and spectral methods. The contour 
levels are the same in the two diagrams* The finite difference calculation 
used a 75 x 50 grid whereas the spectral result was obtained with 
a 32 X 16 grid, with a CFL number of 0*50 and solution smoothing using the 
exponential cut-off applied every 80 time-steps. The major difference between 
the results is that the spectral calculation does not have as deep a pressure 
minimum as the finite difference result* 

Supersonic Flow Past a Circular Cylinder 

The classical problem of a blunt body such as a circular cylinder in a 
supersonic stream has been an ideal test problem for numerical methods since 
it provides a relatively simple well-posed transonic problem with nontrivial 
initial and boundary conditions* The present spectral method obtains the 
steady state solution as the time asymptotic solution of the unsteady Euler 


(33a) 

(33b) 
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equations which are written in the cylindrical polar coordinate (r,0) 

system* The physical domain of interest consists of the known body 

r = r (0), the unknown shock location r = r (0,t), the axis of symmetry (at 
D S 

the front stagnation streamline 0 = tt) and the outflow boundary 

0 = IT - For the purpose of shock fitting, the coordinate 

transformation 


r - 

r (0,t) - r, (0) 
s b 


(34a) 



max 


(34b) 


is introduced so that the shock wave and the body are coordinate lines in the 
transformed domain* The transformed equations of motion, in the notation of 
the shock interaction problems, are 


where 


Qt + ® Qx ^ Qy ^ 


B 


U 

(a^/Y)X^ 

(a^/Y)(l/r)Xg 


YX^ (Y/r)XQ 

U 0 

0 U 


0 


0 0 


0 

0 

0 

u 


(35) 


(36) 


C 


V 

(a^/Y)Y^ 

(a^/Y)(l/r)YQ 


YY^ (Y/r)Yg 0 

V 0 0 

0 V 0 

0 


(37) 


0 


0 


V 
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and 


with 


R 


- lY^ 




U=X^+uX +-X- 
t r r 0 


V 

V = - Y. 
r 0 


(38) 


(39) 


The flow field variables are expanded in double Chebyshev series, and the 
solution technique is the same as for the previous problem* 

The shock boundary r = rg(0,t) (i.e., X = 1) is computed using the 

Rankine-Hugoniot jump conditions and the compatibility equation along the 
outgoing characteristic from the high pressure side of the shock. On the body 
r = (i*e., X = 0) , the normal component of velocity u is zero. The 

limiting angle chosen so that the outflow boundary Y = 1 is 

supersonic, and hence no boundary conditions need be imposed. 

At the symmetry line 0 = tt (or Y = 0) the tangential velocity 

component v is set to zero. The variables P, S, and u (as well as the 
shock velocity) satisfy the condition that their derivatives with respect to 
Y are zero there. This is enforced at each stage of the predictor-corrector 
time discretization (see Eqs. (21) and (22)) by simply using the value zero 
and not the standard Chebyshev spectral Y-derivative values for P, S, and 
u at Y = 0. 

The filtering employed in the calculations reported below was solution 
smoothing every 50 time-steps using the quartic taper (Eq. (12) of Ref. 1) 
with 0^ = 2 tt/ 3. After each filtering step the boundary conditions were 

applied; u was set to zero on the body and v to zero on the symmetry line; 
moreover, the Neumann boundary conditions at Y = 


0 were enforced by 
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transforming to wavenumber space, adjusting the very highest Chebyshev 
coefficient as needed, and then transforming back to physical space. 

Several calculations have been performed for the flow of an intially 
uniform stream past a circular cylinder. The limiting angle was 80®, 

the collocation grid was 9x9, the CFL number was 0.20 and 2000 time-steps 
were taken. Results for the Mach 4 case are illustrated in Fig. 6. Note that 
the essential features of the flow are evident even on this very coarse 
grid. (indeed, it is the small number of data points which is responsible for 
the jagged appearance of the contour lines.) Similar results have been 

obtained for the Mach 2 and Mach 6 situations. Table III presents a 
comparison of the computed values of the stagnation pressure with the 
theoretical results. Since the numerical computations have converged to 
only 3 or 4 digits after 2000 time-steps, the performance of the spectral 

discretization may be even better than implied by Table III. We re-emphasize 
the fact that there remains the clear need for effective means of surmounting 
the severe explicit time-step restriction which besets current Chebyshev 

spectral methods. 

Results for a more challenging flow are shown in Fig. 7. Finite 

difference results (on a 20 x 30 grid) are given in Fig. 8 for comparison. 
The linearly-sheared stream produces a recirculating region. The 9x9 
spectral grid is still capable of resolving this essential feature. 
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Table III. Comparison of Stagnation Pressures for Uniform Flow 
Past a Circular Cylinder 


Mach 

Number 

Calculated 

Pressure 

Theoretical 

Pressure 

Percent 

Relative 

Error 

2 

5.651 

5.6408 

0.18 

4 

21.072 

21 .0750 

-0.014 

6 

46.846 

46.8109 

0.075 


VI. Conclusions 

Our results demonstrate that spectral shock-fitting methods for 
compressible flows are viable techniques. The quantitative comparison for the 
shock/acoustic wave problem shows the superior performance of the spectral 
method. Similar performance is observed on the shock/vortex and blunt body 
problems when the spectral results are compared with finite difference results 
obtained on a much finer grid. 

Our experience shows that before the full potential of spectral methods 
is realized several aspects must be improved. First, filtering techniques for 
both Fourier and Chebyshev methods need to be refined. For filtering in 
Chebyshev methods the problem of conservation and boundary conditions must be 
resolved. Finally, for non-periodic problems the collocation grid 
distribution imposes a severe restriction on the explicit time-stepping used 
throughout this paper. Implicit time-stepping, on the other hand, involves 
expensive inversion of full matrices. There is a clear need to develop 
efficient acceleration techniques. 
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Fig. 1 
Fig. 2 

Fig. 3 
Fig. 4 

Fig. 5 

Fig. 6 
Fig. 7 
Fig. 8 


Figure Captions 


• Typical shock-fitted time dependent flow model in the physical plane. 

• Spectral pressure distribution for a Mach 1.3 Karman vortex street 
showing sensitivity to inflow boundary conditions. 

• Post-shock dependence of the pressure response to a pressure wave 
incident at 10® to a Mach 3 shock. The solid line is the linear 
theory prediction. The circles are the spectral solution. 

• Dependence on incident angle of the pressure response to a 0.1% 
amplitude pressure wave incident on a Mach 3 shock. The solid line 
is the linear theory result. Circles are spectral solutions, squares 
are finite difference solutions. 

• Pressure contours for spectral ( SP) and finite difference (FD) 
calculations of the shock/vortex interaction problem. The spectral 
solution used a 32 x 16 grid and the finite difference used 
a 75 X 50 grid. 

. Spectral solution on a 9 x 9 grid for a circular cylinder in a Mach 4 
uniform stream. 

. Spectral solution on a 9 x 9 grid for a circular cylinder in a 
linearly sheared stream. 

. Finite difference solution on a 20 x 30 grid for a circular cylinder 
in a linearly sheared stream. 
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